Lattice methods for algebraic modular forms on 
classical groups 



Matthew Greenberg and John Voight 



c/3 



Abstract We use Kneser's neighbor method and isometry testing for lattices due 
to Plesken and Souveigner to compute systems of Hecke eigenvalues associated to 
definite forms of classical reductive algebraic groups. 



1 Introduction 

Let Q(x) = Q(x\ ,.. . ,x„) £ 1[x\ , . . . ,x„] be an even positive definite integral quadratic 
form in n variables with discriminant N. A subject of extensive classical study, con- 
tinuing today, concerns the number of representations of an integer by the quadratic 
form Q. To do so, we form the corresponding generating series, called the theta 
series of Q: 

Q (q)= Z q Q[x) eZ[[q]]. 

xeZ n 

By letting q = e 2mz for z in the upper half -plane Jff, we obtain a holomorphic func- 
tion 6 : Jf? — > C; owing to its symmetric description, this function is a classical 
modular form of weight re/2 and level 4N. For example, in this way one can study 
the representations of an integer as the sum of squares via Eisenstein series for small 
even values of n. 

Conversely, theta series can be used to understand spaces of classical modular 
forms. This method goes by the name Brandt matrices as it goes back to early 
work of Brandt and Eichler [17, 18] (the basis problem). From the start, Brandt 
matrices were used to computationally study spaces of modular forms, and explicit 
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algorithms were exhibited by Pizer [42], Hijikata, Pizer, and Shemanske [25], and 
Kohel [36]. In this approach, a basis for S2(N) is obtained by linear combinations 
of theta series associated to (right) ideals in a quaternion order of discriminant N; 
the Brandt matrices which represent the action of the Hecke operators are obtained 
via the combinatorial data encoded in the coefficients of theta series. These meth- 
ods have also been extended to Hilbert modular forms over totally real fields, by 
Socrates and Whitehouse [49], Dembele [8], and Dembele and Donnelly [10]. 

The connection between such arithmetically-defined counting functions and 
modular forms is one piece of the Langlands philosophy, which predicts deep con- 
nections between automorphic forms in different guises via their Galois representa- 
tions. In this article, we consider algorithms for computing systems of Hecke eigen- 
values in the more general setting of algebraic modular forms, as introduced by 
Gross [23]. Let G be a linear algebraic group defined over Q, a closed algebraic 
subgroup of the algebraic group GL„ . (For simplicity now we work over Q, but in 
the body we work with a group G defined over a number field F; to reduce to this 
case, one may just take the restriction of scalars.) Let G(Z) = G(Q) nGL„(Z) be 
the group of integral points of G. 

Suppose that G is connected as an algebraic variety and reductive, so that its 
maximal connected unipotent normal subgroup is trivial (a technical condition im- 
portant for the theory). Let Goo = G(M) denote the real points of G. Then Goo is a 
real Lie group with finitely many connected components. 

Now we make an important assumption that allows us to compute via arithmetic 
and lattice methods: we suppose that Goo is compact. For example, we may take G 
to be a special orthogonal group, those transformations of determinant 1 preserving 
a positive definite quadratic form over a totally real field, or a unitary group, those 
preserving a definite Hermitian form relative to a CM extension of number fields. 
Under this hypothesis, Gross [23] showed that automorphic forms arise without 
analytic hypotheses and so are called algebraic modular forms. 

Let © = Q ®z Z be the finite adeles of Q. Let K be a compact open subgroup of 
G = G(Q) (a choice of level), let G = G(Q), and let 

Y = G\G/K. 

The set Y is finite. Let W be an irreducible (finite-dimensional) representation of G. 
Then the space of modular forms for G of weight W and level K is 

M(W,K) = {f:d/K^W\ f(yg) = yf(g) for all y G G}. 

Such a function / E M(W,K) is determined by its values on the finite set Y; indeed, 
if W is the trivial representation, then modular forms are simply functions on Y. The 
space M(W, K) is equipped with an action of Hecke operators for each double coset 
KpK with p E G; these operators form a ring under convolution, called the Hecke 
algebra. 

Algebraic modular forms in the guise of Brandt matrices and theta series of 
quaternary quadratic forms, mentioned above, correspond to the case where G = 
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PGLi (B) — B x /F x where B is a definite quaternion algebra over a totally real field 
F. The first more general algorithmic consideration of algebraic modular forms was 
undertaken by Lanksy and Pollack [37], who computed with the group G = PGSp 4 
and the exceptional group G = G2 over Q. Cunningham and Dembele [7] later 
computed Siegel modular forms over totally real fields using algebraic modular 
forms, and Loeffler [39] has performed computations with the unitary group U(2) 
relative to the imaginary quadratic extension Q(\/— 11)/Q and 1/(3) relative to 
Q(V— 7)/Q. In this paper, we consider the case where the group G arises from a 
definite special orthogonal or unitary group. Our main idea is the use lattice meth- 
ods, making these computations efficient. This connection is undoubtedly known to 
the experts, and our small contribution is make it explicit and discuss the relevant 
computational aspects. We conjecture that, assuming an appropriate analogue of the 
Ramanujan-Petersson conjecture, lattice methods will run in polynomial time in the 
output size. (This is known to be true for Brandt matrices, by work of Kirschmer 
and the second author [31].) 

To illustrate our method as we began, let Q be a positive definite quadratic form 
in d variables over a totally real field F, and let G = SO(Q) be the special orthogonal 
group of Q over F. (To work instead with unitary groups, we simply work with a 
Hermitian form instead.) Then G is a connected reductive group with = G(F <8>q 
R) compact. Let A be a Z f -lattice in F d . Then the stabilizer K C G of A — A ® z Q 
is an open compact subgroup and the set Y = G\G/K is in natural bijection with 
the finite set of equivalence classes of lattices in the genus of A, the set of lattices 
which are locally equivalent to A . 

The enumeration of representatives of the genus of a lattice has been studied in 
great detail; we use Kneser's neighbor method [34]. (See the beginning of Section 5 
for further references to the use of this method.) Let p C Zp be a nonzero prime ideal 
with residue class field F p . We say that two Zp -lattices A,LT C F" are p-neighbors 
if we have pA , pTI c A n JT and 

dim Fp A / (A fl IT) = dim Fp LT/(A n II) = 1 . 

The p-neighbors of A are easy to construct, are locally equivalent to A , and by strong 
approximation, every class in the genus is represented by a p-neighbor for some p. 
In fact, by the theory of elementary divisors, the Hecke operators are also obtained 
as a summation over p-neighbors. Therefore the algorithmic theory of lattices is 
armed and ready for application to computing automorphic forms. 

The main workhorse in using p-neighbors in this way is an algorithm for isometry 
testing between lattices (orthogonal, Hermitian, or otherwise preserving a quadratic 
form). For this, we rely on the algorithm of Plesken and Souvignier [43], which 
matches up short vectors and uses other tricks to rule out isometry as early as pos- 
sible. This algorithm was implemented in Magma [2] by Souvignier, with further 
refinements to the code contributed by Steel, Nebe, and others. 

These methods also apply to compact forms of symplectic groups; see Chisholm 
[3]. We anticipate that these methods can be generalized to a wider class of reductive 
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groups, and believe that such an investigation would prove valuable for explicit 
investigations in the Langlands program. 

The outline of this paper is as follows. In Section 2, we give basic terminology 
and notation for algebraic modular forms. In section 3, we review orthogonal and 
unitary groups and their Hecke theory. In section 4 we discuss elementary divi- 
sors in preparation for section 5, where we give an exposition of Kneser's neighbor 
method and translate Hecke theory to the lattice setting. In section 6, we present the 
algorithm, and we conclude in section 7 with some explicit examples. 



2 Algebraic modular forms 

In this first section, we define algebraic modular forms; a reference is the original 
work of Gross [23]. 



Algebraic modular forms 

Let F be a totally real number field and let 

F x , = F(g) Q R = M.^. 

Let Q = Q ® x Z be the finite adeles of Q, let F = F <E) @ Q be the ring of finite adeles 
of F. 

Let G be a connected, reductive algebraic group over F. We make the important 
and nontrivial assumption that the Lie group Goo = G(Foo) is compact. Let G = G(F) 
andG = G(F). 

Let p : G — > W be an irreducible (finite-dimensional) representation of G defined 
over a number field E. 

Definition 2.1. The space of algebraic modular forms for G of weight W is 

f ~ /is locally constant and 

M(G,W) = \f:G^W /(7?) = 7/(?)forallreGand?e a }■ 

We will often abbreviate M(W) = M(G,W). 

Each / € M(W) is constant on the cosets of a compact open subgroup K C G, so 
M(W) is the direct limit of the spaces 

M(W,K) = (f.G^W f / (7 ^= 7/ fl p 1. (1) 
v ; \ foi all y e G, g <E G, u £ K J 
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of modular forms of level K. We will consider these smaller spaces, so let K C G be 
an open compact subgroup. When W = E is the trivial representation, M(W, K) is 
simply the space of /s-valued functions on the space Y — G\G/K. 

Proposition 2.2 ([23, Proposition 4.3]). The set Y = G\G/K is finite. 

Let h = #Y. Writing 

h 

G=\jGxiK, (2) 
i=i 

it follows from the definition that any / e M(W,K) is completely determined by the 
elements /(x/) with i = l,...,h. Let 

r i = Gnx i Kxr 1 . 

The (arithmetic) group f/, as a discrete subgroup of the compact group Goo, is finite 
[23, Proposition 1.4]. 

Lemma 2.3. The map 

h 

M(W,K)^®H°(ri,W) 
1=1 

/H. (f(xi ),..., f{x h )) 

is an isomorphism of F -vector spaces, where 

H°(r i: W) = {v e W : yv = vfor all y G I]}. 

In particular, from Lemma 2.3 we see that M(W,K) is finite-dimensional as an E- 
vector space. 



Hecke operators 

The space M(W,K) comes equipped with the action of Hecke operators, defined as 
follows. Let J^(G,K) = J$f(K) be the space of locally constant, compactly sup- 
ported, /iT-bi-invariant functions on G. Then M^iG^K) is a ring under convolution, 
called the Hecke algebra, and is generated by the characteristic functions T(p) of 
double cosets KpK for pGG. Given such a characteristic function T(p), decompose 
the double coset KpK into a disjoint union of right cosets 

KpK = \JpjK (3) 
and define the action of T(p) on / £ M(W, K) by 
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(T(p)f)(g)=Zf{gPj). (4) 

i 

This action is well-defined (independent of the choice of representative p and repre- 
sentatives pj) by the right ^T-invariance of /. Finally, a straightforward calculation 
shows that the map in Lemma 2.3 is Hecke equivariant. 



Level 

There is a natural map which relates modular forms of higher level to those of lower 
level by modifying the coefficient module, as follows [12, §8]. Suppose that K' < K 
is a finite index subgroup. Decomposing as in (2), we obtain a bijection 

h h 
G\G/K' = [J G\(GxiK)/K ' ^ □ rjXKi/Kl 
i=i i'=l 

G(yxiu)K' i y f^S^ 1 )^ 

for y £ G and u£K. This yields 

h 

M(W,K') ^ //°(r?,Hom(^/f;,W)) = 0H°(I-,Coind| W). 

i=i 

Via the obvious bijection 

Ki/K'^K/K', (5) 
letting W = Coindf , W we can also write 

M(W,K' )*±®H°(r i ,W i ) (6) 
1=1 

where W, is the representation W with action twisted by the identification (5). More- 
over, writing K = (Kp)p in terms of its local components, for any Hecke operator 
T(p) such that 

if p <£K' p thentfp = K' p 

(noting that p£K' p for all but finitely many primes p), the same definition (4) applies 
and by our hypothesis we have a simultaneous double coset decomposition 

K'pK' = |J pjk' and KpK = \_\ pjK. 

i i 
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Now, comparing (6) to the result of Lemma 2.3, we see in both cases that modular 
forms admit a uniform description as /t-tuples of ^-invariant maps. For this reason, 
a special role in our treatment will be played by maximal open compact subgroups. 



Automorphic representations 

As it forms one of the core motivations of our work, we conclude this section by 
briefly describing the relationship between the spaces M(W) of modular forms and 
automorphic representations of G. Suppose that W is defined over F (cf. Gross [23, 
§3]). Since Goo is compact, by averaging there exists a symmetric, positive-definite, 
Goo-invariant bilinear form 

(,) : Woo x Woo — >Foo. 
where W» = W ®f Foo. Then we have a linear map 

f :M(W) — ^ Hom Goo (W„o, L 2 (G\(Gx Goo), Foo)) 

by 

n/)(v)(s,S~) = <P(s~)v>/(?)} 

for / g M(W), v g W, and g G x Goo. The Hecke algebra Jt?(K) acts on the 

representation space 

Hom Goo (Wo,L 2 (G\(Gx Goo),/Vo)) 

via its standard action on L 2 by convolution. 
Now, for a nonzero v g W»o, define 

%:M(W)^L 2 {G\(GxGoo),F o) 
/->«P(/)(v). 

(In practice, it is often convenient to take v to be a highest weight vector.) 

Proposition 2.4. The map W r is Jif (K)-equivariant and induces a bijection between 
irreducible J4?(K)-submodules of M(W ,K) and automorphic representations % of 
G(Af ) such that 

(i) 7t(K) has a nonzero fixed vector, and 

(ii) 7Too is isomorphic to poo. 

In particular, an J^(K) eigenvector / g M(W,K) gives rise to an automorphic 
representation. Since automorphic representations are of such fundamental impor- 
tance, explicit methods to decompose M(W,K) into its Hecke eigenspaces are of 
significant interest. 
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3 Hermitian forms, classical groups, and lattices 



Having set up the general theory in the previous section, we now specialize to the 
case of orthogonal and unitary groups. In this section, we introduce these classical 
groups; basic references are Borel [1] and Humphreys [26]. 



Classical groups 

Let F be a field with charF ^ 2 and let L be a commutative etale F-algebra equipped 
with an involution _ : L — > L such that F is the fixed field of L under _ . Then there 
are exactly three possibilities for L: 

1 . L — F and _ is the identity; 

2. L is a quadratic field extension of F and _ is the nontrivial element of Gal(L/F); 
or 

2'. L~FxFand (b,a) = (a,b) for all (a,b) eFxF. 

(As etale algebras, cases 2 and 2') look the same, but we will have recourse to single 
out the split case.) 

Let V be a finite-dimensional vector space over L. Let 

<p : V x V — > L 

be a Hermitian form relative to L/F, so that: 

(i) <p(x+y,z) = q>(x,z) + <p(y,z) for aHx,y,z e V; 

(ii) (p(ax,y) = a(p(x,y) for all x,y E V and a E L; and 

(iii) <p(y,x) = <p(x,y) for all x,y £ V. 

Further suppose that <p is nondegenerate, so (p(x,V) = {0} for x e V implies x = 0. 
For example, the standard nondegenerate Hermitian form on V = L" is 

n 

<p(*,y) = (7) 
i=i 

Let A be the (linear) algebraic group of automorphisms of (V, <p) over F: that is 
to say, for a commutative F-algebra D, we have 

A(D) = Aut D! g, FL (V F ® F D, <p). 

(Note the tensor product is over F, so in particular we consider V as an F-vector 
space and write V F .) More explicitly, we have 

A(F) = Aut L (V,<p) = {T g GL(V) : <p(Tx,Ty) = f(x,y)}. 
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Since <p is nondegenerate, for every linear map T : V — > V, there is a unique linear 
map T* :V ->• V such that 

(p(Tx,y) = (p(x,T*y) 

for all x,y E V. It follows that 

A(F) = {Te GL(V) : TT* = 1} 

where the * depends on (p. 

The group A is reductive but is not connected in case 1 . Let G be the connected 
component of the identity in A. 

In each of the three cases, we have the following description of G < A. 

1. If L = F, then <p is a symmetric bilinear form over F and G = SO (<p) < 0(<p) = A 
are the special orthogonal and orthogonal group of the form <p. 

2. If L is a quadratic field extension of F, then q> is a Hermitian form with respect 
to L/F and G = U((p) = A is the unitary group associated to <p. 

2'. If L = F x F, then actually we obtain a general linear group. Indeed, let e\ = 
(1,0) and <?2 = (0,1) be an F -basis of idempotents of L. Then V\ = e\V and 
V% = ejV are vector spaces over F, and the map T i-> gives an isomorphism 
of G = AontoGL(Vi). 

Remark 3.1. It would also be profitable to consider other groups of symmetries of 
(V, <p), for example, the spin group Spin (<p) in case 1 and the special unitary group 
SU(<p) in case 2. We have simply made one such choice for the purposes of this 
article. 

Remark 3.2. To obtain symplectic or skew-Hermitian forms, we would work instead 
with signed Hermitian forms above. 

Remark 3.3. We have phrased the above in terms of Hermitian forms, but one 
could instead work with their associated quadratic forms Q : V — > L defined by 
Q(y) = (p(v,v). In characteristic 2, working with quadratic forms has some advan- 
tages, but in any case we will be working in situations where the two perspectives 
are equivalent. 



Integral structure 

Suppose now that F is a number field with ring of integers TLp. By a prime of F we 
mean a nonzero prime ideal of Zp. 

Let (V, <p) and G < A be as above. Since our goal is the calculation of algebraic 
modular forms, we insist that Goo = G(Foo) = G(F ®q R) be compact, which rules 
out the case 2' (that L = F x F) and requires that F be totally real. 
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Let 7L L be the ring of integers of L. Let A C V be a lattice in V, a projective Z^- 
module with rank equal to the dimension of V. Suppose further that A is integral, 
so <p(A,A) C Z/,. Define the t/MoZ lattice by 

A* = {ieV : <p(A,x) CZ L }. 

We say A is unimodular if A # = A. 

To a lattice JTCV we associate the the lattice 

n = n®i L z L <zv = v® L t 

with flp = 77 ®z L ^L,p ; we have 77 g)^ Z£ )P = A ®£ L Z/,, p for all but finitely primes 
p. Conversely, given a lattice (77 p ) p C V with flp = A (g>z L Zi jP for all but finitely 
many p, we obtain a lattice 

77 = {jc G V : x G 77 p for all p}. 

(In fact, one can take this intersection over all localizations in V, not completions, 
but we do not want to confuse notation.) These associations are mutually inverse to 
one another {weak approximation), so we write 77 = (JT p )p unambiguously. 
Let 

K = {geG:gA=A} (8) 
be the stabilizer of A in G. Then K is an open compact subgroup of G. Further, let 

r = {geG:gA=A} 

be the stabilizer of A in G. Then the group F is finite, since it is a discrete subgroup 
of the compact group Goo [23, Proposition 1.4]. 

Remark 3.4. In fact, by work of Gan and Yu [21, Proposition 3.7], there is a unique 
smooth linear algebraic group A over with generic fiber A such that for any 
commutative Zp -algebra D we have 

A(D) = Aut D0% Zi (A ® Zl D, <p). 

As we will not make use of this, we do not pursue integral models of A any further 
here. 

We now consider the extent to which a lattice is determined by all of its localiza- 
tions in this way: this extent is measured by the genus, which is in turn is given by 
a double coset as in Section 2, as follows. 

Definition 3.5. Let A and fl be lattices in V. We say A and fl are {G-)equivalent 
(or isometric) if there exists y G G such that yA = fl. We say A and fl are locally 
equivalent (or locally isometric) if there exists g G G such that gA = 77. The set of 
all lattices locally equivalent to A is called the genus of A and is denoted gen(A). 
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For any g = (g p ) p £ G, we have 

& = risp A p ; 

p 

since gpAp = A p for all but finitely many p, by weak approximation, there is a 
unique lattice LT C V such that LT = gA . By definition, LT £ gen(A ) and every lattice 
LT £ gen(A) arises in this way. Thus, the rule 

(g 1 A)^n = gA 

gives an action of G on gen(A). The stabilizer of A under this action is by definition 
K, therefore the mapping 

G/£^gen(A) 
gK^gA 

is a bijection of G-sets. The set G\gen(A) of isometry classes of lattices in gen(A) 
is therefore in bijection with the double-coset space (class set) 

Y = G\G/K. 

By Proposition 2.2 (or a direct argument, e.g. Iyanaga [27, 6.4] for the Hermitian 
case), the genus gen(A) is the union of finitely many equivalence classes called the 
class number of A, denoted h = h(A). 

In this way, we have shown that an algebraic modular form / £ M(W,K) can be 
viewed as a function on gen(A). Translating the results of Section 1 in this context, 
if Ai , . . . , Af, are representatives for the equivalence classes in gen(A ), then a map 
/ : gen(A) — s- W is determined by the finite set /(Aj), . . . ,/(A/,) of elements of 
W. The problem of enumerating this system of representatives for Y becomes the 
problem of enumerating representatives for the equivalence classes in gen(A), a 
problem which we will turn to in Section 5 after some preliminary discussion of 
elementary divisors in Section 4. 



4 Elementary divisors 

In this section, we give the basic setup between elementary divisors and Hecke 
operators, providing a link to the neighbor method in the lattice setting. The results 
are standard. We work in the local case. 

Let F be a local field of mixed characteristic with ring of integers Z/r , let <p be a 
Hermitian form on V relative to L/F, and let Zj, be the integral closure of in L 
with uniformizer P. 

Suppose that G is split, and let T s C G be a maximal split torus in G. Let W s be 
the Weyl group of (G s , T s ). Let 
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X, {T s ) = Uom(T s ,G m ) and X* {T,) = Hom{G m ,T s ) 
be the groups of characters and cocharacters of T s , respectively. 
Theorem 4.1. The Cartan decomposition holds: 

G= \J KX{P)K. 

Xex*(T s )/w s 

There is a standard method for producing a fundamental domain for the acton of 
W s on X*(T S ), allowing for a more explicit statement of the Cartan decomposition. 
Let <P + C X* (T s ) be a set of positive roots and let 

F+ = {Xe X*(T S ) : A (a) > Ofor all a e <£+}. 

Proposition 4.2. Y + is a fundamental domain for the action ofW s on X*(T S ). There- 
fore, we have 

G= \J KX(P)K. 

Agf+ 

We now proceed to analyze this decomposition in Proposition 4.2 explicitly in 
our situation. We suppose that A is unimodular (in our methods, we will consider 
the completions at primes not dividing the discriminant), and we consider two cases. 

We first consider the split case (2') with L = F x F. Let 

V\=e ] _V c ^F n and Ai = eiA^Z£. 

Recall that T n- T\v l identifies G = G(F) with GL(Vi). Already having described 
V\ and A\ in terms of coordinates we have 

G = GL d (F) and K = GL d (Z F ). 

Let T < G be the subgroup of diagonal matrices, consisting of elements t = 
diag(?i ,...,t n ). For i = l,...,n, define the cocharacter X, ■ : F x — > T by 

A,(a) =diag(l,...,l,a,l,...,l), 

where a occurs in the 2-th component. Then 

Y+ = {*?-K':ri<-<r n }. 

Proposition 4.3 (Elementary divisors; split case). Let A and TI be unimodular 
lattices in V with L^FxF. Then there is a basis 

&\ ) • ■ • 3 



of A and integers 



r\<-<r» 
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such that 

(P/PY^e h ...,(F/PY"e„ 

is a basis of LT. Moreover, the sequence r\ < ■ ■ ■ < r n is uniquely determined by A 
and n. 

Now we consider the more difficult cases (1) and (2), which we can consider 
uniformly. We have that either L = F or the maximal ideal of TLp is inert or ramified 
in Z/,. Let v = v(<p) be the Witt index of (V, <p), the dimension of a maximal isotropic 
subspace. Then v = dim 7] <n/2 and V admits a basis of the form 

■ ,e v ,gi, ■ ■ ■ ,gn-2v,fl,- ■ ■ ,fv 

such that 

(p{e h ej) = (p(fijj) = and <p(e,,//) = <5, 7 . 
In this basis, the matrix <p is 

M<P)= ^ (9(gi,8j))tj 

where I is the V x v identity matrix. The set of matrices of the form 
^diag(/i,...,/ v ) 

I 

diag(f I ,...,F v )" 1 

constitute a maximal split torus in G. Considering A, as a cocharacter of GL V (F) as 
above, define 

Mi' :F X ^T 

(Ha) 
a M> jU,(fl) = / 

V A; (op 1 

With these choices, we have 

Y+ = {ii[ l ---ll r f-n<--<r v }. 

Proposition 4.4 (Elementary divisors; nonsplit case). Let A and II be unimodular 
lattices in V with L ^ F x F. Then there is a basis 

ei,---,e v ,gi,...,g„-2v>fi,---,fv 

of A and integers 

n < ■■■ <r v 

such that 
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P n e h ...,P rv e v ,g h ...,gj,p-^f l ,...,P-^f v 

is a basis of 77. Moreover, the sequence r\ < ■ ■ ■ < r v is uniquely determined by A 
and 77. 



5 Neighbors, lattice enumeration, and Hecke operators 

In this section, we describe the enumeration of representatives for equivalence 
classes in the genus of a Hermitian lattice. We develop the theory of neighbors with 
an eye to computing Hecke operators in the next section. 

The original idea of neighbors is due to Kneser [34], who wished to enumerate 
the genus of a (positive definite) quadratic form over Z. Schulze-Pillot [47] im- 
plemented Kneser's method as an algorithm to compute the genus of ternary and 
quaternary quadratic forms over Z, and Scharlau and Hemkemeier [46] developed 
effective methods to push this into higher rank. For Hermitian forms, Iyanaga [28] 
used Kneser's method to compute the class numbers of unimodular positive defi- 
nite Hermitian forms over Z[z] of dimensions < 7; later Hoffmann [24] pursued the 
method more systematically with results for imaginary quadratic fields of discrimi- 
nants d = — 3 to d = —20 and Schiemann [45] extended these computations further 
for imaginary quadratic fields (as far as d — —455). For further reference on lattices, 
see also O'Meara [41, Chapter VIII], Knus [32], Shimura [48], and Scharlau [44]. 



Neighbors and invariant factors 

Let F be a number field with ring of integers Z/r. Let Z/, be a field containing F with 
[L : F] < 2, ring of integers Z/,, and involution _ with fixed field F. In particular, we 
allow the case L~F and Z^ = 7hp . Let (p be a Hermitian form on V relative to L/F, 
and let A cV be an integral lattice. 

If II C V is another lattice, then there exists a basis e\,...,e n for V and fractional 
ideals 2ti , . . . , 2l n and *Bi , . . . , QS n of Zl such that 

A =2liei© •••©2l„e„ 

and 

n = Q3iei©---© $ B„e„ 
(a direct sum, not necessarily an orthogonally direct sum) satisfying 

Q5 1 /2li2--0 5B B /2l„. 

The sequence Q3i/2ti , . . . ,*B„/2l„ is uniquely determined and called the invariant 
factors of fl relative to A . Note that 77 C A if and only if the invariant factors are 
integral ideals of Zl. 
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Define the fractional ideal 

d(A,n) = f\%/% 

1=1 



and let 0(A) = 0(A # ,A), where A # D A is the dual lattice of A. Then in fact 5 (A) = 
0(A), so 0(A) arises from an ideal over %p, which we also denote 0(A) and call the 
discriminant of A . In particular, A is unimodular if and only if (A ) = Zjr , and more 
generaly A p = A ®2 F Z/r p is unimodular whenever p is a prime of F with p { 0(A ). 

Definition 5.1. Let *p be a prime of L and let k G Z with < A: < «. An integral 
lattice L7 C V is a ^-neighbor of A if LT has k invariant factors *p and i.e., 
invariant factors 




k k 

if k<n/2 and 

77 — 2/£— W «— jfe 

if£>n/2andq3^<p. 

We require that ?p 7^ ?P if & > n/2 because the maximal isotropic subspaces in 
this case have dimension < n/2, by Proposition 4.4. It follows from a comparison 
of invariant factors that 77 is a ^-neighbor of A if and only if 

n/(AnU) = (Z L /«P)* andA/(Ann) = (Zl/W- 
A ^-neighbor 77 of A has the same discriminant 0(A) = 0(77). 

Remark 5.2. One may also define A^-neighbors for N a finitely generated torsion 
Zz,-module. 



Neighbors and isotropic subspaces 

Let *p be prime of L above p and let q = Then q = p or q = p 2 , where p is 
the prime below Suppose that p \ 0(A). Let X C A be a finitely generated Z/,- 
submodule. We say thatX is isotropic modulo q if 

<p(x,;y) S q for all x,y e X. 

Define the dual of X to be 

x # = {yey :? (jfj)czi}. 

Then 

¥X # = {yeV:(p(X,y)cy} 
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and A n tyX* C A C V is a lattice. 

Proposition 5.3. Lef ICAte isotropic modulo q. 77zen 

A(«p,x) = qr 1 x+(Anpf # ) 

is a «P* -neighbor of A, where k = AimX/^X. Moreover, A(<$,X) = A($$,X') if 
andonlyifX/y5X=X'/yx'CA/ < $A. 

Proof. The integrality of A (^P,X) is easy to verify using the fact that (p(x,y) <G 
for all x,y e X and <p(x,y) £ for all x el andy e A n^X # . 

First, we prove a claim: A nA(?p,X) = A n^X # . The inclusion Q) is clear. For 
the reverse, suppose y G A HA soy = v+w with v G ^"'X and w G A n<}3X # ; 

then v = y — h> G A and 

9(x,v) c (p(x^- l x) = (p(x,x)y- 1 cyyy- 1 =y 

so v G A PI tyX # and thus y = v + wGAn *}3X # as well. This proves the claim. 
Now, choose a Zz,/<}3-basis x\,...,Xk forX /Cpx C A /*PA. Consider the map 

?(-,*)= A ->(Z L /V)* 

y^(<p(y,x,))+^. 

Since p is coprime to (A ), the Hermitian form <p is nondegenerate modulo since 
X is totally isotropic, it follows that (p(-,X) is surjective. Since (p(y,x) = (p(x,y) for 
all x,y G V, we have that ker<p(-,X) = A n <]3X # . Therefore, by the claim, we have 

A /(A nA(qj,x)) = A/(A r#x # ) = (Z L /^) A \ (9) 

Next, we have 

qp-'xnA =x. 

Therefore, 

A(qj,x)/(AnA(qj,x)) = (qr'x+AnqjxVOAnpf*) 

sqr^/(qrl*n(Anpr*)) (io) 
= «}r 1 x/x^x/q3x 

and X/qJX ^ {1 L /*9) k - Together with (9), we conclude that A(q3,X) is a 
neighbor. 

For the final statement, if X/<}3X = X'/tyX' then it is clear that A(*P,X) = 
A(yi,X') as this construction only depends on X modulo the converse follows 
from the fact that the identification in (10) is canonical: if A PP,X) = A (*$,X') then 

x/qsx=x'/qix'. 

Proposition 5.4. Let n be a ty k -neighbor of A. Suppose that q = is coprime to 
d (A). Then there exists X C A isotropic modulo *P with dimX/*pX = k such that 

n=A(y,x). 
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Proof. Let X be the Z^-submodule of <}37I generated by a set of representatives for 
<pn modulo qp(nn A). Then X is finitely generated and (p(X,X) C = q by the 
integrality of 77. Since 77 is a ^-neighbor, we have 

n/{Ann)^(z L /%) k 

so «pn UnllCA, showing that XCAandX /q3X = (Z L by nondegener- 
acy. 

Next, we prove that 77 C A (^P,X). If y € IT, then by the integrality of 77, 

C 9(fpn,y) = y<p{n,y) C <p. 

Therefore, A n 77 C A n *#X* . But 

a n n c <p _1 x + a n n c 17 

and 

+ (a nn))/(A nn) = qr'x/tqr'xn (a nu)) 
= x/(xn*p(Ann)) ^x/<px 

by construction; since n/(AnU) = (Z L /qj)* andX/^X = (Z L /^) k , we conclude 
s p-'X + (Ann) =77. Thus 

17 = y l x + (a n 77) c qr 1 + (a n^x # ) = A(q3,x) 

as claimed. 

But now since both A(^,X) and 71 are ^-neighbors of A, they have the 
same invariant factors relative to A, so the containment 77 C A(^3,X) implies 
77=A(q3,X). 

Putting together Propositions 5.3 and 5.4, we obtain the following corollary. 

Corollary 5.5. The map X i— > A(^P,X) gives a bijection between isotropic sub- 
spaces X <Z A modulo *p tvif/j dimX/Cpx = A: anc/ ?)3 -neighbors of A. 

From this corollary, we see that by taking a flag inside an isotropic subspace X 
with dimX /tyX = k, every ^-neighbor 71 can be obtained as a sequence 

Aj =A(«p,Xi),A 2 =Ai(?p,X 2 ) J ...,ir = At-i(VA-i) 

of ^-neighbors. However, not all such fc-iterated neighbors are -neighbors: A is 
itself a ^-neighbor of any of its ^-neighbors, for example. 



Neighbors, the genus, and strong approximation 



The ^-neighbors can also be understood very explicitly when *p is odd (i.e., *p \ 2). 
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Let X CAbe isotropic modulo *p with dimX /tyX = k. We revisit the elementary 
divisor theory of Section 4. There is a Z^p-basis xi, . . . ,x n for A p such thatxi,. . . ,x& 
is a basis forX p such that a basis for A (*}2,X)p is 



{P/P)x 1 ,...,(P/P)x k ,x k+u ...,x„ 



if <P^q}and is 



Px\ , . . . ^Px/t^Xk+i , . . . ^X n — k ,P Xn—k+l , . . . ,P X n , 

if fp = fp, where Pisa uniformizer of *p. In the first case, since q = PP has 
qq = 1, this diagonal change of basis is an isometry and thus A($i,X)p = A p ; a 
direct calculation in the latter case shows again that it is an isometry. 

Since the invariant factors of a ^-neighbor are supported over p, we have proven 
the following lemma. 

Lemma 5.6. Let Yl be a ^-neighbor of A with p below *P and p \ 5(A). Suppose 
that p is odd. Then FI belongs to the genus of A. 

Now we form the graph of ^-neighbors: the vertices consist of a set of equiva- 
lence classes of lattices in the genus of A, and for each vertex II we draw a directed 
edge to the equivalence class of each ^-neighbor of fl. This graph is K -regular, 
where K is the number of isotropic subspaces of A p modulo *p of dimension k — 
since all lattices in the genus are isomorphic. If, for example, (p is the standard form 
(so is totally split) and *p ^ then this number is simply the cardinality of the 
Grassmanian Gi(n 1 k)(¥^) of subspaces of dimension k in a space of dimension n, 
and we have the formula 



(q k_ l){q k_ q) ... {q k_ q k-iy 



#GrM)(F ? ) = ^ % & ' (ID 



To conclude, we show that in fact the entire genus can be obtained via iterated 
CP-neighbors; this is equivalent to the assertion that the graph of ^-neighbors is 
connected. 

First, we need the following important result, a consequence of strong approxi- 
mation. For the orthogonal case L = F, see Eichler [16], Kneser [34], or O'Meara 
[41, §104]; for the unitary case L/F, see Shimura [48, Theorem 5.24, 5.27] (and 
Schiemann [45, Theorem 2.10]); and for a further perspectives, see the survey by 
Kneser [35]. 

We say that a lattice A is nice at the ramified primes if for all q ramified in L/F, 
the lattice A q splits a one-dimensional sublattice. If n is odd or L = F or A is even 
unimodular, this condition holds. 

Let Cl(Zi) be the class group of Zi and let Cl(Zi)^ ' be the set of those classes 
that have a representative 21 with 21 = 21. Note C1(Z L ) = C1(Z L )<"> if L = F. 

Theorem 5.7 (Strong approximation). Suppose that 
(i) L = F, « > 3, and 0(A) is squarefree, or 
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(ii) [L : F] = 2, n > 2, and A is nice at the ramified primes. 

Let S be a nonempty set of primes ofL coprime to d(A) that represent all elements 
in Cl(Zt)/Cl(Zt)<~>. 

Then every lattice in gen(A) is equivalent to a lattice LT with LT q = A q for all 
primes q below a prime ^ S. 

Proof (sketch). The hypotheses n > 3 and « > 2 are necessary, as they imply that 
the corresponding spin or special unitary group is simply connected; they also fur- 
ther imply that for all primes p \ J) (A) below a prime *}3 £ S, the form <p p on V p is 
isotropic, so G p is not compact. Since the set S is nonempty, strong approximation 
then implies that every lattice in the spin genus or special genus of A is equivalent 
to a lattice LT as in the statement of the theorem. Finally, the hypothesis that 9(A) 
is squarefree in the orthogonal case implies that the genus of A contains only one 
spinor genus (see O'Meara [41, §102] or Kneser [33]); and the difference between 
the special genus and the genus of A is measured by the group Cl(Zz,)/Cl(Z/,)' ^ 
by work of Shimura when A is nice at the ramified primes. 

Remark 5.8. These are not the minimal set of hypotheses in which strong approxi- 
mation holds, but they will suffice for our purposes; see the references above for a 
more comprehensive treatment. 

We then have the following corollary; see also Kneser [34, §2], Iyanaga [28, 
2.8-2.1 1], and Hoffmann [24, Theorem 4.7]. 

Corollary 5.9. Under the hypotheses of Theorem 5.7, every lattice in gen(A) can 
be obtained as a sequence of ^-neighbors for *}} £ S. 

Proof. Let LT £ gen(A). By strong approximation, we may assume that LTq = Aq 
for all ^ S. 

First, suppose that LTc> = Aq for all ^ Then *}3'"LT C A for some m £ Z>q. 
We proceed by induction on m. If m = 0, then LT C A and since S(LT) = 5(A) we 
have LT = A . 

Suppose m > 0, and choose m minimal so that ^'"LT C A. Let X be the Z/,- 
submodule of ^'"LT generated by a set of representatives for *P"'LT modulo *p m LT n 
*PA. Then X C ^'"LT C A. Now consider again the proof of Proposition 5.4. We 
see that X is isotropic (in fact, (p(X,X) £ q" 1 ). Form the neighbor A($$,X). Then 
A can be obtained from a sequence of ^-neighbors of A. 

Now we have 

$pm-i n = + ( A n $pm-ij7) c qr ! x + (A nyx # ) = A{ty,X) 

since 

<p{X,y, m - l n) C < p(sp m LT,q3 m - 1 LT) C q3q n, ->(TT,TT) C <p. 

Therefore, by induction, <P" LT can be obtained by a repeated ^-neighbor of 
A and we are done by transitivity. 

In the general case, we simply repeat this argument for each prime *p in Z/,. 
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We now connect the theory of neighbors to Hecke operators via elementary divisors 
and the Cartan decomposition as in the previous section. Specifically, we compute 
the action of Jif(G(F p ),K p ) on M(W,K) for primes p \ 0(A). In this case, the cor- 
responding lattice is unimodular. 

As should now be evident from this description in terms of maximal isotropic 
subspaces, the Hecke operator acts on a lattice by a summation over its neighbors. 
We record this in the following theorem. 

Theorem 5.10. Let p€G correspond to the sequence 



in Proposition 4.4 or 4.3. Write KpK = |J PjK. Then for anyx^z G, the set of lattices 



Proof. Each 17, is indeed a ^-neighbor, as is visible by looking at the correspond- 
ing invariant factors. Since the -neighbors are in bijection with maximal isotropic 
subspaces and so are the cosets pj, these sets are in bijection. 



6 Algorithmic details 

Having discussed the theory in the previous sections, we now present our algorithm 
for using lattices to compute algebraic modular forms. 



General case 

We first give a general formulation for algebraic groups: this general blueprint can 
be followed in other situations (including symplectic groups, exceptional groups, 
etc.). We compute the space M(W,K) of algebraic modular forms of weight W and 
level K on a group G. To begin with, we must decide upon a way to represent in bits 
the group G, the open compact subgroup K, and the G-representation W so we can 
work explicitly with these objects. Then, to compute the space M(W,K) as a module 
for the Hecke operators, we carry out the following tasks: 

1. Compute representatives x/K (i — l,...,h) for G\G/K, as in (2), compute f; = 
GnxiKxf 1 , and initialize 



0< ••■ <0< 1 < ■■■ < 1 




Tlj = xpjA — xpjA n V 



is in bijection with the set of ffi -neighbors ofxA. 
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h 

H = ®H°(T i ,W). 

(=1 

Choose a basis of (characteristic) functions / of H. 
2. Determine a set of Hecke operators T(p) that generate J4?(K), as in Section 4. 
For each such T(p): 

a. Decompose the double coset KpK into a union of right cosets pfK, as in (3); 

b. For each % and pj, find yn G G and j* so that 

XipjK = YijXfK. 

c. Return the matrix of T(p) acting on H via the formula 

(7W)(£)=L^/m(%) 
j 

for each / in the basis of H. 

In step 2c, since each function f m is a characteristic function, we are simply 
recording for each occurrence of j* an element of G. 

We now turn to each of the pieces of this general formulation in our case. 



Representation in bits 

We follow the usual algorithmic conventions for number fields [5]. A Hermitian 
form (V, <p) for L/F is represented by its Gram matrix. We represent a Z/r-lattice 
A C V by a pseudobasis over Zp , writing 

A = 2lixi © • • • © 2l„x„ 

with xi,... ,x n E V linearly independent elements and 21, C L fractional Z^-ideals 
[6]. The open compact subgroup K is the stabilizer of A by (8) so no further speci- 
fication is required. 

The irreducible, finite dimensional representations of G are given by highest 
weight representations. The theory is explained e.g. by Fulton and Harris [19], and 
in Magma there is a construction of these representations [13, 4], based on the LiE 
system [38]. 



22 



Matthew Greenberg and John Voight 



Step 1: Enumerating the set of representatives 

We enumerate a set of representatives x,K for G\G/K using the results of Sections 
4 and 5. For this, we will use Corollary 5.9, and so we must assume the hypotheses 
of Theorem 5.7, namely: 

(i) L = F, n > 3, and 0(A) is squarefree, or 

(ii) [L : F] = 2, n > 2, and A is nice at the ramified primes. 

Next, according to Corollary 5.9 we compute a nonempty set of primes S of L co- 
prime to 25(A) that represent all elements in Cl(Z/,)/Cl(Z/,)' '. By the Chebotarev 
density theorem, we may assume that each prime is split in L/F if L ^ F . There 
are standard techniques for computing the class group due to Buchmann (see Co- 
hen [6, Algorithm 6.5.9] for further detail). We compute the action of the involution 
_ on C1(Z/,) directly and then compute the subgroup Cl(Zi,)' ) fixed by _ and the 
corresponding quotient using linear algebra over Z. 

Next, we traverse the graph of CP-neighbors for each *}5 £ S. To do this, we per- 
form the following tasks: 

a. Compute a basis for Am as in Propositions 4.3 and 4.4 according as L = F or 
L + F. 

b. Compute the one-dimensional isotropic subspaces modulo *}3 in terms of the basis 
e,- for the maximal isotropic subspace. 

c. For each such subspace X, compute the ^-neighbor A($5,X) = ?fi~ l X + ^X # 
using linear algebra. 

d. Test each neighbor A for isometry against the list of lattices already com- 
puted. For each new lattice A', repeat and return to step a with A' in place of 
A. 

Since the genus is finite, this algorithm will terminate after finitely many steps. 

Remark 6.1. One can also use the exact mass formula of Gan and Yu [21] and Gan, 
Hanke, and Yu [20] as a stopping criterion, or instead as a way to verify the correct- 
ness of the output. 

Further comments on each of these steps. 

First, in steps la-lb we compute a basis. When [L: F]= 2, this is carried out as 
in Section 5 via the splitting L p = F p x F p . When L = F, we use standard methods 
including diagonalization of the quadratic form: see e.g. work of the second author 
[51] and the references therein, including an algorithm for the normalized form 
of a quadratic form over a dyadic field, which at present we exclude. From the 
diagonalization, we can read off the maximal isotropic subspace, and this can be 
computed by working not over the completion but over If /p e for a large e. Next, in 
step lc we compute the neighbors. This is linear algebra. Step Id, isometry testing, 
is an important piece in its own right, which we discuss in the next subsection; as a 
consequence of this discussion, we will also compute 17 = Aut(A,). From this, the 
computation of a basis for H = 0f = i H°(Tl,W) is straightforward. 
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To test for isometry, we rely on standard algorithms for quadratic Q-spaces and Z- 
lattices even when computing relative to a totally real base field F or a CM extension 
L/F. Let a\,... ,a,j be a Z-basis for Zi with a\ = 1, and let jci, . . . ,x n be a basis of 
V. Then 

{<M,}/ ! ,/ 

j=l,...,n 

is a Q-basis of V. Define Q-bilinear pairings 

(Pi-.VxV — >Q by <Pi[x,y) = tr L/Q (p(a iX ,y). 

Since a i = 1 and <p is a definite Hermitian form on V over L, <pi is a positive definite, 
symmetric, bilinear form on V over Q. In other words, (V, <pi ) is a quadratic Q- 
space. The L-space (V, <p) can be explicitly recovered from (V, <pi), together with 
the extra data (fa, . . . , (pj by linear algebra. Note that the forms (p2 , . . . , q>d are in 
general neither symmetric nor positive definite. 

Lemma 6.2. Let A and FI be lattices in V. A *L-linear map f : A — > LT is an 1l- 
linear isometry if and only if each <p,- is invariant under f. 

Using Lemma 6.2, we reduce the problem of testing if two Hermitian lattices 
over are isometric to a problem of testing if two lattices over Z are isometric 
in a way which preserves each <p,-. For this, we rely on the algorithm of Plesken 
and Souvignier [43], which matches up short vectors and uses other tricks to rule 
out isometry as early as possible, and has been implemented in Magma [2] by 
Souvignier, with further refinements by Steel, Nebe, and others. 

For each representative lattice found, we compute a set of invariants to quickly 
rule out isometry whenever possible. These invariants include things like the sizes 
of the automorphism groups, the first few terms in the theta series, and invariants of 
sublattices (e.g. those generated by short vectors). 

Remark 6.3. An essential speed up in the case of Brandt modules is given by 
Dembele and Donnelly [10] (see also Kirschmer and the second author [31, Algo- 
rithm 6.3]). To decide if two right ideals I, J in a quaternion order are isomorphic, 
one first considers the colon ideal (7 : J)l = { a € B : aJ C 1} to reduce the problem 
to show that a single right ideal is principal; then one scales the positive definite 
quadratic form over Q by an explicit factor to reduce the problem to a single short- 
est vector calculation. It would be very interesting to find an analogue of this trick 
in this context as well. 



Step 2: Hecke operators 

Essentially all of the work to compute Hecke operators has already been set up in 
enumerating the genus in Step 1 . The determination of the Hecke operators follows 
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from Sections 4 and 5, and their explicit realization is the same as in Step la. We 
work with those Hecke operators supported at a single prime. In Step 2a, from The- 
orem 5.10, the double coset decomposition is the same as set of ^-neighbors, which 
we compute as in Step 1 . In Step 2b, we compute the isometry yu using isometry 
testing as in the previous subsection: we quickly rule out invalid candidates until the 
correct one is found, and find the corresponding isometry. Finally, in Step 2c, we 
collect the results by explicit computations in the weight representation. 



7 Examples 

In this section, we illustrate our methods by presenting the results of some explicit 
computations for groups of the form G = U L //r(3), relative to a CM extension L/F, 
where L has degree 2, 4, or 6. We compute the Hecke operator at unramified, degree 
one primes p of L corresponding to summing over p-neighbours of a lattice. 

Remark 7.1. We made several checks to ensure the correctness of our programs. 
First, we checked that matrices of Hecke operators for p and q with p ^ q commuted. 
(They did.) Additionally, a known instance of Langlands functoriality implies that 
forms on U{\) x t/(l) x t/(l) transfer to 1/(3). Checking that resulting endoscopic 
forms occur in the appropriate spaces [39, §§4.2, 4.6] also provided a useful test of 
our implementation. (It passed.) 



Example 7.2. U L/F {3), L = Q(y/^7), F = Q, weights (0,0,0) and (3,3,0): 
Table 1 Computation of T p i onM(Q) for unramified, degree one p C Z L with 2 < N(p) < 200. 



N(p) 


2 


11 


23 


29 


37 


43 


53 


67 


71 


79 


107 


time (s) 


0.02 


0.07 


0.18 


0.28 


0.42 


0.57 


0.82 


1.35 


1.46 


1.86 


3.73 


a p 


7 


133 


553 


871 


1407 


1893 


2863 


4557 


5113 


6321 


11557 


bp 


-1 


5 


41 


-25 


-1 


101 


47 


-51 


185 


-15 


293 






N(p) 


109 


113 


127 


137 


149 


151 


163 


179 


191 


193 


197 


time (s) 


4.18 


4.59 


5.85 


7.08 


8.56 


9.04 


10.88 


13.78 


16.92 


17.22 


17.29 


a p 


11991 


12883 


16257 


18907 


22351 


22953 


26733 


32221 


36673 


37443 


39007 


bp 


215 


-109 


129 


-37 


335 


425 


237 


-163 


-127 


131 


479 



Here, we extend aspects of the calculation in the principal example of [39]. In 
this case, the class number of the principal genus of rank 3 Hermitian lattices for 
L/Q is 2, with classes represented by the standard lattice A\ = I? L and the lattice 
A2 C L? with basis 



(1-0,0,0), (1,1,0), i(-3 + co, -1 + co,-l + co). (co = i(l + V = 7)) 
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The lattices A\ and A2 are 2-neighbours: 

(0A2 C Ai, 0)A\ c Ai. 

(Representatives for the ideal classes in the principal genus were computed in the 
first place by constructing the 2-neighbour graph A\ .) It follows that the space of 
M(Q) of algebraic modular forms for U L / F (3) with trivial coefficients is simply the 
2-dimensional space of Q-valued function on {[Ai], [A2]}. We obtain two distinct 
systems a p and b p of Hecke eigenvalues occurring in M(Q). (See Table 1.) We point 
out observations of Loeffler considering the nature of the corresponding algebraic 
modular forms: First, observe that the system a p is "Eisenstein", in the sense that 
the eigenvalues of T p .1 is the degree of the Hecke operator T p y. 

a p = N{p) 2 +N{p) + l. 

Equivalently, the corresponding algebraic modular form is the lift from 1/(1) x 
t/(l)xt/(l)of Xtrin x Xtriv x Xtriv ■ The algebraic modular form with system of eigen- 
values b p is also a lift form i/(l) x 17(1) x U(l): 

b p =p 2 + pp + p 2 . 

Here, we abuse notation and write p for either of its generators, and note that the 
expression on the right in the above is independent of this choice. 

We now consider analogous computations involving forms of higher weight. The 
space M(\^ 3 0) the associated to the above data has dimension 4, while the represen- 
tation space V3 3 itself has dimension 64. Loeffler [39] showed thatM(V3 3 0) splits 
as the direct sum of two 2-dimensional, Hecke-stable subspaces not diagonalizable 
over(Q>(\/^7): 

M(V XX0 ) =Wi®W 2 . 

One of these spaces arises as the lift involving a 2-dimensional Galois conjugacy 
class of classical eigenforms in Sg{I\ (7)) via a lifting from the endoscopic subgroup 
U(l) x U (2) of 1/(3). The other corresponds to a Galois conjugacy class of nonen- 
doscopic forms, whose associated £-adic Galois representations p : Gz, — s- GL2(Qc) 
are irreducible. 

We consider the corresponding modular space M(Va o/¥-]) of mod (V —7)- 
modular forms of weight (3,3,0). We computed the Hecke operators T p ,\ for un- 
ramified, degree one primes p C Zl of norm at most 100. In Table 2, we present the 
corresponding run-times. 

Simultaneously diagonalizing the matrices of these Hecke operators we obtain 
two mod (\/— 7) systems of eigenvalues that we write a p and b p . We choose this 
notation because those systems are the reductions modulo 7 of the corresponding 
trivial weight systems a p and b p from earlier. We have an explicit modulo 7 congru- 
ence between a nonendoscopic form in weight (3,3,0) and an endoscopic form in 
weight (0,0,0). Thus, the modulo 7 Galois representation associated to the system 
b p is reducible. 
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Table 2 Computation of T p \ mod (V— 7) for degree one, unramified pcZi with 2 < N(p) < 100. 





2 


11 


23 


29 


37 


43 


53 


67 


71 


79 


Steps 2a,b 


0.11 


2.49 


12.37 


26.52 


60.08 


128.29 


265.47 


595.90 


984.19 


1561.67 


Step 2c 


0.41 


26.77 


123.82 


208.38 


328.95 


450.68 


686.25 


1431.75 


1414.790 


1774.52 


d p 











3 





3 








3 





b p 


6 


5 


6 


3 


6 


3 


5 


5 


3 


6 



Example 7.3. U L/F {3), F = Q( vT3), L = F^-13-2VUj , weight (0,0,0): 

In this example, the class number of the principal genus is 9, as is the dimension 
of the corresponding space of automorphic forms with trivial weight. We computed 



Table 3 Computation of Tn\ for degree one, unramified pcZ[ with 2 < N(p) < 250. 



N(p) 


29 


53 


61 


79 


107 


113 


131 


139 


157 


191 


211 


time (s) 


15.15 


51.73 


70.35 


123.21 


216.82 


242.20 


339.50 


378.81 


486.22 


727.81 


943.89 


dp 


4 


11 


12 


5 


5 


3 


6 


10 


1 


11 


2 



the matrices of the Hecke operators acting on the (Q-vector space M(Q) for unrami- 
fied, degree one pcZt with 2 < N(p) < 250. There appears to be a 1 -dimensional 
"Eisenstein" subspace on which T p i acts via degT p i = N(p) 2 + N(p) + 1. The 8- 
dimensional complement of this line decomposes into Q-irreducible subspaces of 
dimensions 2, 2, and 4. In all of these computations, the level subgroup is the stabi- 
lizer of the standard lattice 1? L C L? . 

The Hecke algebra does acts nonsemisimply on the space M(Fi3) of modulo 
13 automorphic forms. It appears that the minimal polynomial of T v i has degree 
6 when N(p) = 1 (mod 13) and degree 7 when N(p) = 3,9 (mod 13). (Other 
residue classes do not occur for norms of degree one primes of F splitting in L.) 
When N(p) = 1 (mod 13), the eigenvalue 3 = JV(p) 2 +N(p) + 1 (mod 13) oc- 
curs with multiplicity 5, while when N(p) = 3,9 (mod 13), the eigenvalue = 
N(p) 2 +N{p) + l (mod 13) occurs with multiplicity 1. Finally, is a 1 -dimensional 
eigenspace in M(Fi3) with eigenvalues a q as in Table 3. 

Example 7.4. U L/F (3), K = Q(C 7 + ^'),L = Q(£ 7 ), weight (0,0,0): 

In this case, the class number of the principal genus is 2, with the classes repre- 
sented by the standard lattice A\ = 1? L and its 29-neighbour with basis 

(C7 + C7-C 7 5 ,0,0), (6,1,0), 

^(-138 - 234^7 - 210C 2 - 303C 3 - 258£ 4 - 1 17C 5 , 
16 + 12C+13C 2 + 20C 3 + HC 4 + 6C 5 , 

16 + 12C+ 13C 2 + 20C 3 + HC 4 + 6C 5 )- 
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Table 4 Timings: computation of T p A mod (■\/ = 7) for p with 2 < N(p) < 100 and p ^ p 



N(p) 


29 


43 


71 


113 


127 


197 


211 


239 


281 


time (s) 


2.16 


2.85 


6.77 


16.43 


21.35 


51.14 


53.58 


73.05 


101.84 


a P 


871 


1893 


5113 


12883 


16257 


39007 


44733 


57361 


79243 


bp 


-25 


101 


185 


-109 


129 


479 


-67 


17 


395 



Automorphism groups of A\ and Ai and their roles in this computation. As 
above, a p =N(p) 2 +N(p) + 1 = degTn i, and the form with system of Hecke eigen- 
values fl p is a lift from U(l) x U(l) x U (I). Also, observe that 

flp = bp = 3 (mod 7), 

implying that the modulo 7 Galois representation attached to the system bp is re- 
ducible. 
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